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I.  INTRODUCTION 

This  report  addresses  Itself  to  the  question  of  what  limits  sound  speed 
fluctuations  in  the  ocean  place  on  the  dimensions  and  integration  times  of 
acoustic  arrays.  We  shall  for  the  present  defer  any  discussion  of  advanced 
processing  techniques,  and  confine  ourselves  to  the  simple  case  of  an  array 
that  measures  incident  Intensity  from  various  angles  by  adding  linearly  varying 
phase  shifts  to  the  hydrophones  in  the  array,  and  by  integrating  the  received 
signal  over  some  period  of  time.  We  investigate  at  what  array  dimensions  and 
at  what  integration  time  coherent  processing  is  destroyed  by  sound  speed 
fluctuations  Induced  by  the  ocean  medium.  The  results  obtained  for  this 
situation  will  serve  as  a convenient  baseline  against  which  to  measure  the 
improvements  that  may  be  obtainable  through  the  use  of  more  sophisticated 
processing  techniques. 


II . COHERENT  PROCESSING  AND  THE  PHASE  STRUCTURE  FUNCTION 

Let  us  assume  that  we  have  a receiver  consisting  of  a vertical  rectangular 
array  of  hydrophones  of  horizontal  and  vertical  lengths  y and  z with  a density 
of  n hydrophones  per  unit  area.  Let  us  further  assume  that  there  is  a source 
located  at  a horizontal  bearing  angle  9,  (relative  to  the  array  normal)  emitting 
an  acoustic  frequency  o^.  Then  in  the  absence  of  perturbations  Induced  by  the 
ocean,  the  receiver  would  pick  up  a signal 


N , , Y Z T Iqy (sine  -sine)  -l(o  -o)t 
I / dy  / dz  / dt  e ° e ° 


(2.1) 


a=l 


0 


0 


0 


if  it  Integrates  for  a time  T at  a frequency  a and  is  looking  in  the  horizontal 

direction  0.  Here  a indexes  the  different  acoustic  rays  along  which  the  signal 

(a) 

can  travel  from  source  to  receiver  and  \p  is  the  complex  jignal  associated 
with  ray  a.  The  acoustic  wave  number  q is  defined  by  q = c/Oq  where  c is  tVie 
sound  speed  at  the  receiver. 
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In  general,  the  ocean  medium  through  which  the  signal  travels  will  disturb 
the  received  signal  so  that  the  signal  associated  with  each  ray  path  a will 

Cq.) 

vary  from  one  hydrophone  to  another,  and  will  vary  with  time.  Thus  ijj  ^ becomes 

replaced  by  ^(x,t),  and  hence,  when  averaged  over  an  ensemble  of  oceans,  the 

received  intensity  from  direction  0 at  frequency  a will  be 

„ Y Y Z Z T T 
(1(6, a))  = n \ / dy  / dy'f  dz  / dz'f  dt  Jdt' 

a a'  0 0 0 0 0 0 

iq(sin0^-alne)(y-y')  -i(»„-o) (t-t ' ) *,(«'),*,  ,,, 

e e (x,t)  ip^  (x',t')> 

(2.2) 

The  fluctuations  induced  by  the  ocean  along  two  different  rays  are  inde- 
pendent. Hence  ip^°‘'  for  a ^ a' . These  averages  are 

-1/2  -1/2 

proport  'nal  to  e e (see  Ref.  1).  When  the  fluctuations 

2 2 
(a)  (a') 

are  lai  ge  $ and  $ are  large  and  terms  with  a a'  contribute  little 


to  the  received  intensity. 

For  terms  with  a - a'  we  have 

,(a)  * 


1) 


{^P'"^'(x,t)  ip' 


(a) 


where  D is  the  phase  structure  function  and  is  the  intensity  average  of 

ray  a.  To  a good  approximation  D is  Independent  of  ray  and  we  shall  assume  this 


henceforth.  Thus 


(1(0,0))  = n“  ^ 1^“^  / dy  / dy'  I dz  I dz'  J dt  j , 


a 


0 


0 


0 


0 


0 


dt' 


iq(sin0Q-sin0) (y-y')  -i(Op-o) (' -t' ) ~h  D(x-x',  t-t') 
e e e 

(2.3) 

and  our  problem  is  reduced  to  the  study  of  the  behavior  of  D. 

Before  proceeding  to  do  this,  however,  let  us  anticipate  what  we  will 


learn  about  D to  give  an  estimate  of  when  coherent  processing  fails. 
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First,  we  note  that  as  either  the  spatial  or  temporal  separations  go  to 
2 2 

infinity,  D ->■  2$  where  $ is  the  parameter  defined  in  Ref.  1,  and  which 

coincides  with  the  mean  square  phase  fluctuation  along  a ray  when  the  fluctuation 
2 

is  small.  Thus  $ sets  the  basic  scale  of  D. 

2 

If  $ is  sufficiently  small,  then  D is  everywhere  small  and  there  is  no 
limit  on  coherent  processing.  To  make  a numerical  estimate  of  when  this  occurs, 
we  note^^  that  for  an  axis  ray 


where  R is  the  propagation  range  and  v is  the  acoustic  frequency.  The  region 

2 

in  R and  v where  $ <1,  where  coherent  processing  is  always  valid,  is  shown 

in  Fig.  1. 

2 

When  R and  v are  such  that  >1,  limits  on  coherent  processing  become 
meaningful.  In  this  case  only  small  spatial  or  temporal  receiver  separations 
are  relevant  and  (as  will  be  shown  later)  D is  essentially  quadratic  in  the 
separations.  Hence,  taking  for  simplicity  the  case  of  a point  receiver  inte- 
grating over  time,  we  may  set 


D(0,t-t’)  = 


-) 


(2.4) 


where  t is  some  characteristic  time.  With  this  form  for  D,  we  may  now  evaluate 

the  integral  in  (2.3)  (for  0=  0^,  a = o^)  and  we  find 

T T - k D(0,t-t') 

<I(T))  = / dt/  dt'  e 

0 0 


T ^ 


T large 
T small 


(2.5) 


A schematic  plot  of  the  integral  as  a function  of  T is  shown  in  Fig.  2.  The 
dotted  line  represents  the  true  behavior,  interpolating  smoothly  between  the  two 
asymptotic  expressions  given  in  eq.(2.5).  A reasonable  definition  of  where 
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coherent  processing  fails  is  where  the  two  asymptotic  forms  meet;  this  occurs 
at  T = xi/2^/0.  At  this  value  of  T,  d(0,T)  = 2ti . Thus  we  shall  use  D = 2tt 

as  our  criterion  for  the  failure  of  coherent  processing. 

III.  GENERAL  FORMULAE  FOR  THE  PHASE  STRUCTURE  FUNCTION 

The  calculation  of  the  phase  structure  function  D can  be  carried  out  as  follows 
ihe  most  general  situation  is  one  in  which  we  are  comparing  two  rays  joining  two 
distinct  sources  to  two  distinct  receivers,  (Fig.  3).  We  shall  label  the  rays 

1 and  2;  ray  1 joins  a source  at  to  a receiver  at  and  ray  2 connects  a 

source  at  x^2  to  a receiver  at  Xj^2*  let  x^  and  X2  be  horizontal  distance  along 
rays  1 and  2 respectively;  then  (x^^.y^^Cx^) , z^(x^))  and  (x2,y2(x2),  Z2(x2>)  denote 
points  on  the  two  rays.  We  therefore  have  , Zj^(Xgj^))  = (Xgj^.ygj^.Zg^) 

and  I “ ^\l*^Rl’^'Rl^  ’ analogous  expressions  for  ray  2. 

The  phase  structure  function  for  these  two  rays  is  defined  by 

D(l,2)  ^ ds  - q / ^ dq)2>  . (3.1) 

ray  1 ray  2 

In  Ref.  1 it  is  shown  that  when  the  two  rays  are  not  far  apart,  and  when  the 
sound-speed  fluctuations  are  dominated  by  internal  waves,  (3.1)  may  be  replaced  by 

D(l,2)  = 2 / dx  / dk  / dm 

0 0 0 ^y 

[1  - cosAt  cos(k  Ay(x)  + k Az(x))]  . (3.2) 

y 2 

where  Ay(x)  and  Az(x)  are  the  horizontal  and  vertical  separations  of  the  two 
rays  at  a horizontal  distance  x along  the  rays,  and  At  = tj^~t2  is  the  time 
separation  of  the  rays,  and  F(k,m;  z(x))  .'.s  the  internal-wave  spectrum.  The 
spatial  separations  of  source  and  receiver  come  in  through  the  boundary  conditions 
on  Ay(x)  and  Az(x): 


> 


4y(0)  = yg^  - yg2  = Vg,  Sy(R)  - - y^^  = \ 


Az(0)  = Zg^  - Zg2  + tan6(0) (Xg^-Xg^)  = ^g  + tan6(0)Xg  , 


Az(R)  = - Zj^2  tane  (R)  (Xj^j^-Xj^2)  = ■*■  tan6(R)Xj^  . 

Since  the  deterministic  sound  speed  c(z)  varies  only  vertically,  we  have  an 
explicit  expression  for  Ay(x),  namely: 


S 


Ay(x)  = Yg  + (Yj^  - Yg)  x/R 

The  quantity  Az(x)  is  defined  as  the  solution  of  the  differential  equation 


(3.3) 


Az(x)  + 1^-^  c(z(x)^  Az(x)  = 0 , 


(3.4) 


dx  c *-dz 

subject  to  the  boundary  conditions  given  above. 

As  an  illustration,  suppose  the  sound  channel  is  parabolic:  c(z)  = 
2-2 

c + 1/2  K (z-z)  , where  z is  the  sound  channel  depth.  Then  the  differential 
equation  for  Az(x)  has  the  solution 


Az(x)  = 


(Z_+tan6(R)X_)slnKx  + (Z  +tan0 (0)X„) sinK(R-x) 
K K D b 


sin  KR 


(3.5) 


More  generally,  the  solution  to  tV»e  differential  equation  for  Az(x)  can 
be  written 

^ dx’ 

Az(x)  = atan9(x)  + 3tan6(x)  J 2 > 

0 (tan6(x)) 

where  the  constants  a and  3 are  to  be  determined  from  the  boundary  conditions 
given  earlier.  For  very  long  ranges,  when  there  are  many  loops  in  the  ray,  the 
integral  over  dx'  becomes  nearly  proportional  to  x.  Thus,  we  have,  approximately 

Az(x)  = tan6(x) 


p.Zg  + tanS(0)Xg. 

A , /^R 

L\  tan0(O)  / \ ^ 

/ V tan0(R)  /U/_ 

(3.6) 
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The  quantities  k and  k in  (3.2)  are  given  by  the  internal-wave  dispersion 

y 2 


relations: 


k = k 

y 


/2 

/O) 

2 

- “l 

2 

0) 

•H 

3 

1 

(3.7) 


k = k 
z 


n(z) 


,/  2 2 
\m  “ 


(3.8) 


2 2 2 2 

where  O)^  = o)  + n 0 . Here  n(z)  is  the  buoyancy  frequency  and  o)  is  the 
inertial  frequency.  The  Garrett-Munk  75  model  of  the  internal-wave  spectrum  is 


F(k.g);  z(x))  ^ g ^^2^ 


' (I  -T 

\j=i  j + 


B n 


i ‘ 
4 


0 


2 2 

/2 

2 

3 

1 

3 

/O) 

-0), 

1 

1 

/ 

i 

3 ' 

k 

/ 2 

2 

0) 

/ 0) 

-“t 

.222 

0)  -0), 


0)  -0)^ 
\ B j “"n^ 


2 2 
0)  -0), 


. /O)  -0).  V 

e (n-0))  0 (u)-w^)  0 ^k  - ^ j 


(3.9) 


This  spectrum  is  normalized  by  the  requirement  that  the  parameter  $ 


is 


= / dx  / dk  / d»  . 


0 0 0 


(3.10) 


- 


'?f.  -T*'-  •/ 


1 Vga 


Here  (Hq)  is  the  mean  square  sound  speed  fluctuation  at  the  surface 
(actually  at  the  bottom  of  the  mixed  layer),  is  a mode  number  cutoff,  q is 
the  (sound  axis  value  of  the)  acoustic  wave  number,  n^  is  the  buoyancy  frequency 
on  the  sound  axis,  n^  is  the  buoyance  frequency  at  the  surface,  co^  is  the 

inertial  frequency  and  B it,  the  vertical  scale  depth.  Canonical  numerical  values 
of  these  are: 


(p  ,)  = 4.9  X 10  ^ 

0 rms 


n^/2Ti 


j*  = 3 


0^/271 


m^/2TT 


B 


1.1  cph 
3 cph 
1/24  cph 
1 km 


The  0 functions  in  (3.9)  limit  the  ranges  of  integration  in  (3.10),  and 
in  particular  exclude  singularities  at  m = 0 and  k = 0. 

Expressions  similar  to  (3.10)  obtain  if  we  use,  for  example,  m and  mode 
number  j as  independent  variables: 

R 


D(l,2)  = 2 / dxl  /dm 

0 j=l  0 K 


[1  - cosmAt  cos  (k  Ay(x)  + k Az(x))] 

y z 


(3.11) 


where  now  we  replace  k and  k bv 

y z ^ 


2 2 

. i/U) 

k = JJL  lZ_ L 

y B n. 


0 


n 


Again  the  normalization  is  such  that 


= r V r j..  f^,j;z(x)) 


C3.12) 


and,  for  the  GM75  spectrum,  we  have 


! t 


F_(m,j;z(x))  ^ _8_  2,  __ 

\c  0 “ 

y “ I 


- 42^.2 

j=l  j +J* 


2 3/22 

q B OJ^  n /o3  -0) 

~2 — I 0 (n-oj)  e (w-ojj^) 
0)  -0) 

JLi 


2 3 

Hq  u 


(3.13) 


IV.  INDIVIDUAL  SMALL  SEPARATIONS  AND  THEIR  CONTRIBUTIONS  TO  D 

We  are  primarily  interested  in  cases  where  the  parameter  is  large, 

corresponding  to  the  saturated  or  partially  saturated  regimes,  for  only  here 

will  there  be  limitations  on  coherence.  For  this  case  we  only  need  evaluate 

the  phase  structure  function  for  small  spatial  and  temporal  separations,  since 

the  contribution  of  large  separations  in  eq.(2.3)  is  heavily  damped  by  the  large 
2 

value  of  <D  . When  the  separations  are  small,  their  effect  in  D is  additive,  so 
we  may  discuss  each  of  them  separately.  We  shall  find  that  the  general  formulae 
for  D,  given  in  the  previous  section,  become  relatively  simple  when  the  separations 
are  small,  and  we  will  end  up  with  a rather  simple  analytic  expression  for  D 
characterized  primarily  by  coherence  lengths  for  temporal  and  vertical  and  hori- 
zontal spatial  separations.  Let  us  outline  how  this  comes  about. 

A.  Temporal  Separations 

The  general  expression  for  D reduces  to 


D(T)  = 2 f dxf  dk  /do)  (1 


0 0 


- coswT) 


(4.1) 


For  small  T we  replace  1 - coswT  by  1/2  w^T^  so  that 


D(T)  = $^T^ 


(4.2) 


8 


where  $ is  defined  In  Ref.  1.  Finally,  we  write  this 


as 


D(I)  . . (i  , 

With  v,n-  coherence  time  t defined  simply  by 


(4.3) 


$ 

T = — 


For  an  axis  ray,  this  time  Is^^ 


(4.4) 


2 2 

T = 2ui^  log  n^/oj^  = 0.45  (hours) t = 1.6  hours  . 

B-  Vertical  Spatial  Separations 

Here  It  Is  most  convenient  to  use  the  w-j  language: 

' hk  I (1  - . 

(Z,  the  length  of  a vertical  receiver,  comes  Into  this  formula  through  the 
boundary  conditions  on  Az(x)).  Here 


(4.5) 


k = il  g(z(x)) 


B n 


0 


(4.6) 


The  value  of  using  m and  J as  variables  Is  that  ..ends  only  on  J , and  In 

fact  Is  simply  linear  In  j . Because  of  this  fact,  and  because  F(»,J;t(x))  la 

proportional  to  + J^),  the  sum  on  J can  be  approximately  evaluated 

analytically,  we  replace  J by  /dj , and  define 

j 


I(“»j*)  - / dj  J k y cos  aj/  / dj  ^ ^ 

J A ^JLA  ^ 


i 


0 


j +j: 


jf 


- 2 Cl  (ajg) 


+ [cosh  aj*(Cl(cx(jQ+lj^))  +Cl(a(jQ-lj^)))] 
+ l[slnh  aj*(Si(«(jQ+lj^))  - Sl(a(jQ-lj^)))] 


(4.7) 


where  Cl  and  SI  are  the  cosine  and  sine  Integral  functions.  We  can  set  1 =1 

® -^0 
here,  to  correspond  to  J ' 

j=l 
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II)  1 ' I ti  ■ii^#i  I ‘^rr*-ir  T 


In  terms  of  tuc  function  I,  the  phase  structure  function  can  now  be  written 

R 


»(Z)  - 2 / ■>>■  iz  / .i-  (i  - I'i  aM 


(4.8) 


A simplified  fr>j.m’'la,  valid  for  small  separations,  is  obtained  by  expanding 
the  function  I for  small  Az : 


D(Z)  - 2 / dx  tn  d»  ffl'ciog  f + V 1 , (4.9) 

0 j y V v/  V ^ 


where  is  an  x-dependent  length  defined  by 


I/log  (i^fl)  "o 


£ (x)  = — jQ££(ji±ll 

V 2tt  n(z(x)) 


(4.10) 


and  where  is  a constant  given  by 


= 2 - ^ 


(4.11) 


log(j;+l) 

where  C - 0.577...  is  Euler's  constant. 

The  expression  (4.9)  for  D can  be  evaluated  analytically  for  a near  axis 
ray.  For  these  rays  the  sound  channel  is  nearly  parabolic,  and  the  rays  are 


therefore  nearly  sinusoidal.  Thus  (cf  eq.(3.5)  with  6 = 0 = Z.  = 0 and  = Z) 

K D o K 


Az(x)  = Z 


sin  Kx 


sin  KR 


(4.12) 


Using  this  formula  in  (4.9)  , and  noting  that  F(k,a);z(x))  is  not  a function  of 
X for  the  axis  ray,  we  find 


D(Z)  = 2$2  (i  (logZ+r)}  , 

Av  ^ 


(4.13) 


where  the  coherence  length  '3^^  for  this  ray  is 
Av  v 

(^„  here  is  £ evaluated  on  the  sound  axis)  pnd  the  constant  T is 

V V \r 


■■v  ■ I + I'v  - I '“8  2 . 


Numerically,  we  find 

= 300  meters, 


and 


r = 1.49. 

V 


Eq.(4.9)  for  D can  also  be  evaluated  analytically  for  a steep  ray,  because 
for  such  rays  the  major  contribution  to  the  integral  comes  from  the  region  near 

the  ray  apex. 

For  a single  apex,  we  can  replace  the  quantity 
2 


2U  ) 

V V 


Az  , V 
Og  T"  + Y„) 


in  the  Integrand  by  its  value  at  the  apex.  The  remaining  integral  then  gives 

2 


just  the  parameter  2$^  evaluated  in  the  apex  approximation.  Thus 

2 


D(l,2)  = $ 


(t)' 


(log  + Y^) 


(4.14) 


where  the  hat  on  Az  and  indicate  that  they  are  to  be  evaluated  at  the  ray 
apex  (i,z).  (Note  that  y^  is  Independent  of  (x,z),  so  we  do  not  have  to  indicate 

any  special  point  at  which  it  is  to  be  evaluated.) 

For  very  long  ranges,  where  a steep  ray  has  many  apexes,  we  ^an  use  (3.7) 
to  calculate  Az.  For  example,  for  a point  source  and  a vertical  receiver  of 

height  Z, 


. _ 0 (x)  . xZ 


where  we  have  replaced  tan0  by  0 since  all  angles  are  small,  and  where  Z^^  is  Z. 


With  this  approximation  we  cannot  simply  evaluate  Az  at  the  apex,  because 


vanishes  there.  Rather  we  must  average  - 6^  over  the  ray  apex.  This  yields 


the  expression 


D(Z)  = $ 


2 1 


0^ 


1 2 
^ 0(R) 


(t)^ 


(log  y-  + Y^) 


(4.15) 
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where  the  average  of  0 over  an  apex  is 


6^  .,m  ^ 

1/  TT^  *^0 


"O  ^ 
n^(z)r 


(4.16) 


where  r is  the  radius  of  curvature  of  the  ray  at  the  apex.  An  obvious  modi- 
fication of  (4.15)  gives  the  expression  for  D for  arbitrary  source  and  receiver 
separations  as  in  eq.(3.7). 

We  may  put  eq.(4.14)  into  a form  parallel  to  eq.(4.12)  by  choosing 

^ = /3  fc  , (4.17) 

‘^v  V -2 

as  the  coherence  length  for  a steep  ray,  and 

r = Y + log  . (4.18) 

V V 

Numerically,  for  the  surface  limited  ray,  = 420m  and  = 3. 

C . Horizontal  Spatial  Separations 
Now  we  have 

R 00 

D(Y)  = 2 / dx  y / do)  - (1  - cos  k Ay(x))  , (4.19) 

0 j=l  0 y ^ 


where 


Ay(x)  = Y x/R 


for  a horizontal  receiver  of  length  Y.  Furthermore, 


k = 

y B 


2 2 
■“l 


n 


0 


Unfortunately,  in  contrast  to  case  B,  k^  is  now  not  a function  of  j alone.  An 
analytic  evaluation  of  (4.19)  can  therefore  not  be  carried  very  far.  It  is 
still  possible  to  carry  out  the  sum  on  j,  obtaining  an  analog  of  (4.8);  however, 
the  0)  integration  is  now  prohibitively  difficult.  As  a result,  simple  expressions 
analogous  to  (4.9)  and  (4.13)  are  not  now  obtainjible. 

12 
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Nevertheless,  as  we  shall  see  In  the  following  section,  a fairly  good 
(30%)  analytic  representation  of  D(Y)  Is 


D(Y)  = 2$^  I 


where  the  horizontal  coherence  length  Is  6.4  km. 

SUMMARY 

Combining  these  separations  we  have  the  following  simple  approximate 
expression  for  U,  valid  for  small  separations. 


(i  . I (f)^  <10.^.  Vl 


with 


T = 1.6  hours 


and  for  near  axis  rays 


= 300  m, 

- 6.4  Km. 

and  =1.49 

while  for  the  surface  limited  ray. 


sfv  ■ ““ 


r - 3. 

V 


V.  HORIZONTAL  ARRAYS  AND  STATIONARY  POINT  SOURCES 

As  an  application  of  the  results  obtained  so  far,  let  us  consider  the  case 

of  a stationary  point  source  and  a horizontal  linear  receiver.  Thus  the  source 

separations  are  all  zero  and  are  zero  as  well.  Then  Az(x)  = 0 and 

Ay(x)  = Yx/R  where  we  write  Y for  Y_.  From  eq.  (3.2)  we  have 

K 

D(Y.T)  = 2 / dx  / dk  / do) 

0 0 0 *^y 

1 2 2 

tl  - coswT  cos  (kJ-2 — 2 ^ ^ > (5.1) 

]j  U)  -U)^ 

where  T = t^-t^.  To  illustrate  what  this  implies  we  shall  specialize  to  the 
axis  ray  z(x)  = 0.  For  this  ray  tj  =<*).,  the  integral  on  dx  can  be  carried 
out , and  we  find 


D(Y,T)  = 2R  / 
0 


dk  / do) 
0 


F(k,co;  0) 
k 

y 


[1 


- costoT 


sin  kY 


kY 


(5.2) 


The  expression  for  D(Y,T)  can  be  easily  evaluated  numerically.  Fig.  4 

2 

shows  the  ratio  D(Y,T)/24>  for  various  values  of  T as  a function  of  Y.  Fig. 5 
shows  the  same  ratio  as  a function  of  T for  various  values  of  Y. 

Figures  4 and  5 also  show  lines  representing  fits  to  the  exact  expression 
of  the  form 


D(Y,T)  _ 1 T.2  .1  , Y ,2 

9 o 2 •' 


24) 


H 


(5.3) 


The  value  t 1.6  hours  (which  we've  seen  in  Section  IV  is  the  value  of  4>/$)  is 
imposed  on  the  fit.  The  value  of  is  selected  to  optimize  the  fit  and  turns 
out  to  be  6.4  km. 
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We  may  obtain  an  analytic  result  for 


with  coherent 


maximum  values  of  T and  Y consistent 


processing  by  using  eq.  (5.3); 


(7)^  + (/-)2  = 2TT 


'H  $ 
For  zero  Integration  time  this 


(5.4) 


yields  a maximum  array  length  of 


MAL  = 3^ 


which  for  a near-axis  ray  gives 


(5.5) 


MAL  = 210  \ MOOO  km 


R 


similarly  a maximum  Intagratlon  tlma  can  ba  obtained  an 


(5.6) 


MIT  = r - ^ 

$ ^ ~ V 


(5.7) 


where 


■t, 


H 


V “ = 4 km/hr 


(5.8) 


For  intermediate  oomblnatlona  of  array  length  and  Integration  time  llmlta 
may  be  obtained  from  Figures  4 and  5. 


For  example,  an  acoustic  frequency  of  100  Hz  . 

y luu  Hz  (i.e.,  an  acoustic  wavelength 


. — j waveiengl 

of  15  m)  at  a range  of  1000  km  yields  = 26  2 Tho  t n/o  2 

> as  ^ Zb. 2.  The  line  D/2$^=  27t/(56.4)  = 


.112  interaeota  the  1 . 0.1  hours  curve  In  Pig.  . at  V-  3.3  hm.  Thus  lor  an 
ntegratlon  time  of  0.1  hours  the  maximum  useful  array  length  Is  3.5  km. 


VI.  MOVING  SOIIRrK.t; 

AS  a final  example  of  the  use  of  the  result  obt.,  red  above,  let  us  discuss 
the  effects  of  moving  sources.  Por  a point  source  moving  „lth  velocity 

''  " ‘■Vy’V’  “=  PAplace  the  source  dimensions  (Xg.Yj.Z  ) by  (V  T,v  T.V  T) 

We  shall  confine  ourselves  to  ths  r-no  n ^ y » 2 • 

to  the  case  of  sources  moving  In  the  horizontal  plane 


and  set  = 0.  Be  may  then  distinguish  two  possibilities. 
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^ Source  mover  tranrverse  to  the  Hue  Joining  nource  to  receiver. 

s case  we  may,  as  In  Section  IV,  limit  ourselves  to  discussing  the 
-1=  ray.  „e  then  have  ha(r,  - 0.  and  ay(r,  = V^T  + (v-v  t)x/R  where,  as  In 

Section  IV,  y Is  the  horlaontal  dimension  of  the  receiver.  Thus.  eq.  (3.2) 
becomes  replaced  by 


D(Y,T)  = 2R  / dk  / dm  .O). 

r\  ^ V 


cosmT 


^sln  kY  - sin  kV  T 

; JL_ 

kY  - kV  T 
y 


For  a point  receiver  (V  » 0)  „e  therefore  find  exactly  the  same  results  as 
those  given  In  Section  V,  but  with  V replaced  by  V^T  so  that  the  maxim™ 
Integration  time  In  this  case  Is  (using  our  analytic  approximation) 


MIT  = 


(6.2) 


Thus  the  velocity  of  the  moving  source  adds  quadratlcally  to  the  Intrinsic 
internal-wave  velocity,  v,  to  reduce  the  coherent  Integration  time.  The  value 
Of  at  Which  source  action  begins  to  affect  coherent  Integration  Is  evidently 
Vy  ■v  V - 4 km/hr.  For  Y ^ 0 evaluation  of  e,.(6.1)  must  be  done  numerically. 

£1— Source  moves  directly  toward  or  away  from  receiver. 

We  again  specialise  to  a point  receiver.  Then  dy(x)  - 0 and  In  calculating 

dz(x)  the  boundary  conditions  are  ds(0)  = ta„e(0)VJ,  ds(R)  = 0.  For  short 

integration  times  I we  can  add  the  direct  effect  of  T and  the  Induced  Is  due 

to  the  source  motion.  The  effect  of  is  Is  given  by  eq.  (4.8). 

Insofar  as  the  analytic  formula  (4.17)  is  vaHH  nno  w 

VH.j./;  IS  valid  one  can  estimate  the  maximum 

integration  time  T for  a point  receiver  ^n  Unra 

in  the  same  way  as  in  case  A,  with  Z 

In  (4.17)  replaced  by  V^T9(0);  T satisfies  the  equation 


9 *2  2 

27T  = $ + X 


V T0(O) 


+ VJ  . 


(6.3) 


I* 


nis 


To  estimate  the  value  of  at  whleh  radial  source  motion  begins  to  cut  into 

coherent  integration  time,  we  may  Ignore  the  log  term.  The  relevant  velocity 
is  then  evidently 


/r  0(o)t 


.75  km/hr 


for  the  surface  limited  ray.  This  velocity  is  considerably  smaller  than  the 
analogous  tangential  velocity,  due  to  the  fact  that  the  shorter  vertical  coherence 
length  («0m  as  compared  to  6.4  km)  more  than  compensates  for  the  Increase  of  the 
critical  velocity  by  tbe  factor  1/6(0)  v 5.  Thus  for  steep  rays  radial  source 

motion  is  a much  more  serious  Inhibition  to  coherent  integration  than  is 
tangential  motion. 

In  conclusion  we  should  emphasize  that  the  limitations  produced  by  source 
motion  discussed  here  are  only  those  associated  with  internal  wave  Induced 
fluctuations  in  the  ocean.  There  are  also  effects  produced  simply  by  source 
motion  in  the  presence  of  many  deterministic  multipaths.  These  have  recently 


been  analyzed  by  H.  Cox^'.  (These  determrnistic  multipath  effects  are  of  course 


not  relevant  for  a stationary  source.) 
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Figure  Captions 


Fig.  1;  The  parameter  $ as  a function  of  propagation  range  and  acoustl*' 
frequency. 


Fig.  2:  Schematic  plot  of  received  intensity  at  a point  receiver  as  a function 

of  integration  time,  showing  the  transition  between  the  region  where 
coherent  processing  applies  (I(T)  T^)  and  where  it  fails  (I(T)  “ T) . 

The  time  Is  the  Intersection  of  the  small  T and  large  T asymptotic 

forms  for  I(T).  We  say  coherent  processing  is  valid  below  this  time. 

Fig.  3:  Sketch  illustrating  ray  1,  ray  2,  the  mean  ray,  and  the  source  and 

receiver  separations  projected  in  the  vertical  plane.  A similar  sketch 
(with  straight  line  rays)  applies  projected  ia  the  horizontal  plane. 

2 

Fig.  4:  Graphs  of  D(Y,T)/2$  computed  numerically,  for  various  integration 

times  T as  a function  of  horizontal  receiver  length  Y,  for  a flat  ray 
and  a stationary  point  source.  The  dotted  lines  show  an  approximate 
fit  of  the  form  D/2$^  = (T/t)^  f ~ with;(^  = 6.4  km  and 

T = 1.6  hrs. 


Fig.  5: 


Same  as  Fig.  4 but  plotted  as  a function  of  T for  various  values  of  Y. 
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